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Abstract 



The Regge Calculus is a powerful method to approximate a continuous manifold by a 
simplicial lattice, keeping the connectivities of the underlying lattice fixed and taking 
the edge lengths as degrees of freedom. The Discrete Regge Model limits the choice 
of the link lengths to a finite number. We examine the phase structure of Standard 
Regge Calculus in four dimensions and compare our Monte Carlo results with those 
of the Z 2 -Regge Model as well as with another formulation of lattice gravity derived 
from group theoretical considerations. Within all of the three models of quantum 
gravity we find an extension of the well-defined phase to negative gravitational 
couplings and a new phase transition. We calculate two-point functions between 
geometrical quantities at the corresponding critical point and estimate the masses 
of the respective interaction particles. A main concern in lattice field theories is the 
existence of a continuum limit which requires the existence of a continuous phase 
transition. The recently conjectured second-order transition of the four- dimensional 
Regge skeleton at negative gravity coupling could be such a candidate. We examine 
this regime with Monte Carlo simulations and critically discuss its behavior. 



0.1 Introduction 



The construction of a non-trivial quantum theory of gravitation represents one of 
the major open problems in theoretical physics. Approaches based on perturba- 
tive methods failed due to the non-renormalizability of the underlying theory in 
dimensions greater than two. One possibility to deal with the quantization problem 
consists in discretization of space-time and quantization via the path integral, in 
close analogy to quantum field theory on a flat background geometry. There are 
two related schools attempting at a lattice theory of quantum gravity, dynamical 
triangulation and Regge quantum gravity. We concentrate on the latter approach 
and refer to [1] for an excellent review of the former framework. Regge calculus is 
the only discretization scheme reproducing the Bianchi identities of classical general 
relativity and is formulated exclusively in terms of the edge lengths of the lattice, 
independent of any coordinate system. Taking advantage of the freedom in how 
to choose the lattice structure as well as the action, we investigate in this paper 
different formulations, all of them based on the original work by Regge [2], namely 
Standard Regge Calculus (SRC) [3, 4], a Group Theoretical Model (GTM) [5], and 
the Z 2 -Regge Model (Z 2 RM) [6, 7]. Different to ordinary lattice gauge theory, quan- 
tum fluctuations are represented by fluctuations in the edge lengths. In particular 
the Monte Carlo method is a convenient tool for evaluating the discrete Euclidean 
functional integral in a non-perturbative way. It was shown that within a certain 
range of the gravitational coupling the path integral converges to make up a well- 
defined phase of lattice quantum gravity. We present an extension of the phase 
diagram of the models in four dimensions mentioned above. Due to the fact that 
a Wick rotation from the Lorentzian to the Euclidean sector of quantum gravity 
is not feasible, in general the sign of the action in the path integral is not fixed 
a priori. Therefore we investigate the continuation of the well-defined phase into 
the region of negative gravitational coupling. It turns out to be terminated by a 
pronounced transition-like cross-over to an ill-defined phase suggesting to look for 
a continuum limit. In particular the Z 2 RM permits accurate investigations of the 
negative coupling region. Further, we calculate two-point functions to probe the 
existence of massless quanta of the gravitational field. A candidate for a realistic 
quantum theory of gravity should reproduce the expected long-range interaction 
behavior observed in nature. 

0.2 Lattice quantum gravity 

Any smooth <i-manifold can be approximated by appropriately gluing together pieces 
of flat space, called rf-simplices, ending up with a simplicial lattice. We take the 
edge lengths as the dynamical degrees of freedom and leave the triangulation of the 
lattice fixed. Adopting the Euclidean path-integral approach we may write down 
the partition function 




(1) 



with the gravitational action / introduced in more detail below. The functional 
integration extends over the squared edge lengths qi of the links / of the simplicial 
lattice. 

One of the problems with (1) is the ambiguity in performing the link-lengths inte- 
gration. Commonly the measure is written as 

D[q]=T[d qiq r 1 Hq) , (2) 

i 

with T a function of the quadratic edge lengths being equal to one if the Euclidean 
triangle inequalities are fulfilled and zero otherwise. The question remains whether 
such a local measure is sufficient at all and how the power a might be chosen. In 
an effort to clarify the role of the measure the conventional definition of diffeomor- 
phisms has been employed in the two-dimensional case, assuming that a piecewise 
linear space, i.e. a Regge surface, is exactly invariant under the action of the full 
diffeomorphism group [8] . After a conformal gauge fixing was performed in the con- 
tinuum formalism with the DeWitt measure, it was shown that the evaluation of the 
non-local Faddeev-Popov determinant by using such a Regge regularization leads to 
the usual Liouville field theory results in the continuum limit. All that is based 
on a description of piecewise linear manifolds with deficit angles, not edge lengths, 
and is mostly taken as an argument that the correct measure of Standard Regge 
Calculus has to be non-local. However, to our knowledge it is not obvious that this 
argument carries over to a discretized Lagrangian, which is formulated in terms of 
fluctuating edge lengths, obeying triangle inequalities, and which is not invariant 
under the diffeomorphism group due to the presence of curvature defects: different 
assignments of edge lengths correspond to different physical geometries, and as a 
consequence there are no gauge degrees of freedom in Standard Regge Calculus, 
apart from special geometries (like flat space) [9, 10]. Furthermore, the generaliza- 
tion of this procedure to higher dimensions and its numerical implementation are 
technically demanding and some of the proposed non-local measures do not agree 
with their continuum counterparts in the weak field limit, which is a necessary con- 
dition for an acceptable discrete measure [11]. But this property is fulfilled for the 
standard simplicial measure (2). 

Working in Euclidean space, i.e. with positive definite metric, the conformal mode 
renders the four-dimensional continuum Einstein-Hilbert action unbounded from 
below. This unpleasant feature persists in the discretized Regge-Einstein action 
but need not necessarily lead to an ill-defined path integral [3]. Indeed, numerical 
simulations reveal the existence of a well-defined phase with finite expectation val- 
ues within a certain range of the bare gravitational coupling. A lattice action for 
gravitation is given by 

I(q) = -PY. R t(A,5 t ) + \J2V s , (3) 

t s 

where the first sum runs over all triangles t with areas A t (q) and the corresponding 
deficit angles S t (q) yield the curvature elements R t . The second term extends over 
the volumes V s (q) of all four-simplices s and allows together with the cosmo logical 
constant A to set an overall scale in the action. Next we define our three different 
models for subsequent numerical treatment. 



0.2.1 Standard Regge Calculus (SRC) 



Here we employ the Regge-Einstein action with R t = A t 5t, set the cosmological 
constant A = 1 and choose the gravitational measure to be uniform, a = 1 [3, 4]. In 
the classical continuum limit the Regge action 2 J2t Rt is equivalent to the Einstein- 
Hilbert action / d^x^fgR if the so-called fatness 4> s of a 4-simplex obeys [12] 

A lower limit / = 10~ 4 restricts the configuration space and thus facilitates numer- 
ical simulations [13]. 



0.2.2 Group Theoretical Model (GTM) 

Constructing the dual of a simplicial lattice, Poincare transformations can be as- 
signed to its links to yield an action in which the sin of the deficit angle enters, 
Rt = A t sm5 t , and that reduces to the Regge action in the small curvature limit [5]. 
While in the classical continuum limit the deficit angles become small, this is in 
general not the case in the path-integral quantization on a finite lattice, where one 
has to sum over all possible discrete configurations. A priori it is hence not clear 
whether both SRC and the GTM result in the same quantum continuum theory. 
For comparison we again use the uniform measure with a = 1, set the cosmological 
constant A = 1, and choose a lower limit on the fatness, / = 10~ 4 . 



0.2.3 Z 2 -Regge Model (Z 2 RM) 

This model was invented in an attempt to reformulate SRC as the partition function 
of a spin system [6, 7]. It is defined by restricting the squared link lengths to take 
on only two values 

qi = 6j(l + eai) , a t E Z 2 . (5) 

By setting bi = 1,2,3,4 for edges, face diagonals, body diagonals and the hyper- 
body diagonal of a hypercube, respectively, the link lengths are allowed to fluctuate 
around their flat-space values. The Euclidean triangle inequalities are automatically 
fulfilled as long as e < e max G M + and therefore T = 1 in any case. Consequently, the 
uniform measure in the quantum-gravity path-integral becomes identical to unity 
for all possible link configurations. The action can be rewritten in terms of com- 
plicated local "spin-spin" interactions. Nevertheless, numerical simulations of the 
Z 2 RM become extremely efficient by implementing look-up tables and a heat-bath 
algorithm. Computations have been performed with the parameter e = 0.0875 and 
the cosmological constant A = because (5) already fixes the average lattice volume. 



0.3 Phase structure 



Monte Carlo simulations with at least 50k sweeps for each value of the coupling 
f3 have been performed on regularly triangulated hypercubic lattices with toroidal 
topology and N = 4 4 vertices. We measured expectation values of the average link 
length (q) = (iVf 1 J2i Qi) (Ni denotes the total number of links) and the average 
curvature 



within the models described above. 

Let us first concentrate on SRC, cf. figures la and 2a. The expectation values of the 
average link length change only slightly, increase near (3~ ~ —0.16 and /3+ « 0.116, 
and become arbitrarily large for (3 < (3~ and (3 > (3+ indicating the emergence of 
an ill-defined phase with spike-like structures. The expectation value of the average 
curvature is negative in the well-defined phase, but | (R) | becomes very large in the 
ill-defined phase where simplices collapse into degenerate configurations with long 
links and small volumes. The coupling of SU(2)-gauge theory to SRC was examined 
in [14] and turns out to have little influence on the gravitational phase diagram. 
The GTM behaves very similar, see figures lb and 2b. Again expectation values are 
finite and well-defined only in a certain coupling region between f3~ ~ —0.115 and 
/3+ ~ 0.14. In contrast to the Regge model the curvature (R) is mostly positive in 
this coupling interval. However, as long as we deal with non-renormalized entities 
it is not clear whether the observed differences in the phase structure are physically 
significant. Although SRC and GTM have the same classical continuum limit, it is 
necessary to demonstrate that they coincide in their critical properties in order to 
have the same quantum continuum limit. 

The Z 2 RM is always well-defined because (5) limits the link lengths a priori, cf. 
figures lc and 2c. The lattice freezes, forming characteristic configurations where 
the systems with continuously varying edge lengths become ill-defined. Such a model 
is of course better suited for investigations of the phase boundaries which here occur 
at j3~ ~ —4.665 and /3+ pa 22.3. While the transition at (3+ is quite clearly of first 
order, the behavior at (3~ is more subtle and will be discussed in some detail below 
in Section 0.5. Different to SRC and the GTM the four kinds of links in the Z 2 RM 
are not equivalent by definition (5). This is reflected in the corresponding average 
link lengths q in a single configuration, see figure 3. Thus we can examine their 
respective influence on the phase transitions. For example the hyperbody diagonal 
shows a considerable discontinuity at f3+ and exactly this diagonal was identified 
from weak-field calculations of Regge gravity to belong to the (five) spurious metric 
degrees of freedom [15]. It would be interesting to clarify whether physical and 
spurious degrees of freedom decouple as predicted from weak-field Regge theory. A 
simulation of the Z 2 RM with the hyperbody diagonal fixed results in unchanged 
order of the transition but in slightly shifted values of the critical couplings. An 
intriguing question remains: to what extent are the (ten) physical degrees of freedom 
responsible for (the order of) the phase transitions? 

A related issue was investigated in [16] using the SRC on general, non-regular tri- 




(6) 



angulations of the four-torus. Also there, even with additional higher-order terms 
in the action, the phase transition is influenced by the local coordination numbers. 



0.4 Correlation functions 



One important feature of a physically relevant theory of quantum gravity should be 
the existence of a massless graviton. The gravitational weak-field propagator can 
be cast into a spin- zero and a spin- two part [17], with the latter written in terms of 
connected curvature-curvature correlations, 



The local operators should be measured at two vertices Vq and Vd separated by the 
geodesic distance d, which we take to be equal to the index distance along the main 
axes of the skeleton. This seems a reasonable approximation in the well-defined 
phase with its small average curvature. In general one expects for (7) at large 
distances the functional form 



A power law with a = 2 and a vanishing effective mass to = would hint at 
Newtonian gravity with massless gravitons. Indeed this was found in [17] for SRC 
with higher-order curvature terms at some positive critical coupling. In contrast to 
that, fast decaying exponentials were reported for (7) and also for volume correla- 
tions in the case of SRC at /3+ [18]. 

We now compare SRC and GTM (on 4 3 x 8 lattices) with Z 2 RM (on 8 4 lattice) 
[19]. The gravitational couplings were chosen close to the negative critical coupling 
f3~ where there is a chance for a continuous phase transition and thus a continuum 
limit. Figure 4 displays our Monte Carlo data for the two-point functions (7). We 
took only distances d > 2 into account because Gn(d = 1) is plagued by lattice 
artifacts due to contact terms. In order to test whether (7) obeys a power law we 
fixed a = 2 and fitted the effective masses to, cf. table 1. For SRC to remains rather 
constant towards the critical coupling whereas in the GTM the mass decreases for 
(5 — > f3~. However for all fits the uncertainties in the mass parameters are large 
allowing even for to = 0. For both the Z 2 RM and the GTM, power-law fits to (8) 
are compatible with vanishing effective mass to. Also for SRC an algebraic decay is 
compatible but, due to the large uncertainties, obviously no definite conclusion can 
be drawn here. 




(7) 
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0.5 Continuum limit at negative gravitational cou- 
pling 

The Discrete Regge model Z 2 RM - like full Regge theory - exhibits two phase 
transitions [20]. One is located at a negative value and the other one at a positive 
value of the bare gravitational coupling. Although earlier work concentrated on 
the latter transition [21], there is no reason for favoring a positive value of (3 over 
a negative one due to the Wick rotation problem mentioned in the Introduction. 
In the vicinity of the transition at positive (3, histograms, e.g. of A t S t [20], clearly 
show a two-peak structure, see figure 5. The two phases coexist and tunnelling 
from one phase to the other and back takes place. The system also exhibits a 
hysteresis; the transition occurs at a larger value of f3 if the simulation is started 
from a configuration in the well-defined phase than it does if the calculation is started 
from a "frozen" configuration. Given all these pieces of evidence, we conjecture that 
the transition at (3 > is of first order. 

To determine the order of the transition at negative f3 we employed in Ref. [20] 
histogram techniques and used the Binder-Challa-Landau (BCL) cumulant criterion 
[22]. The BCL cumulant is defined as B L := 1 — 3 |^ 2 | 2 , with E being the action of 
the system under consideration. It was evaluated for A t 5 t on different lattice sizes 
with L = 3 to 10 vertices per direction, simulating the system at several values of 
the bare coupling f3 with high statistics (using always e = 0.0875 and A = 0). For 
the BCL cumulant a trend towards 2/3 was observed and all histograms showed a 
clear one-peak structure, cf. Ref. [20]. In our recent Monte Carlo simulations [23] we 
typically generated 200 000 — 500 000 iterations, and recorded for every run the time 
series of the energy density e = E/N and the magnetization density m = J2i^i/N , 
where N = L 4 is the lattice size. To obtain results for the various observables O 
at values of the bare gravitational coupling (3 in an interval around the simulation 
point Po, the reweighting method [24] was applied. By this means we can compute 
the specific heat, 

C(/3) = /3 2 iV ((e 2 )-( e ) 2 ), (9) 
and the (finite lattice) susceptibility, 

Xifl) = N ((m 2 ) - (H> 2 ) , (10) 

in a certain /5-range around the simulation point. Figure 6 shows the finite-size 
scaling (FSS) of the maxima of the specific heat C max and the susceptibility Xmax, 

Table 1: Effective masses m of the curvature two-point function Gr((T) for gravita- 
tional couplings close to (3~ . 





SRC 






GTM 




Z 2 RM 


(3 -0.155 
m 2(9) 


-0.157 
3(24) 


-0.159 
3(17) 


-0.110 
1.0(21) 


-0.111 
0.9(8) 


-0.112 
0.5(7) 


-4.668 
0.5(9) 



respectively. While the behavior of C max could still be explained as critical scaling 
with a negative exponent a, the flattening of Xmax as L increases is quite unusual for 
a second-order transition and should rather be taken as an indication for a cross-over 
regime. 

Another feature of the system can be seen in figure 7 depicting the critical gravita- 
tional coupling, as determined from the specific-heat maxima C max and the suscep- 
tibility maxima Xmax- In the case of a second-order transition the extrapolations of 
all pseudo-transition points lead to one infinite-volume critical value. This seems to 
be violated in the four-dimensional Discrete Regge model and thus again favors the 
interpretation as a cross-over phenomenon over a true, thermodynamically defined 
phase transition. 

0.6 Summary and conclusion 

We explored the complete phase structure of three different formulations of lattice 
quantum gravity in four dimensions including the region of negative gravitational 
couplings [20]. The qualitative resemblance of the results from simulations of the 
Z 2 -Regge Model to those with continuously varying edge lengths is particularly 
remarkable. We computed two-point functions close to the critical bare gravitational 
coupling in the negative coupling regime. In the case of the Group Theoretical 
Model as well as for the Z 2 -Regge Model we found some evidence for long-range 
correlations corresponding to massless gravitons. Eventually, if there is universality 
between SRC, GTM, and Z 2 RM, the same results for physical quantities are to be 
expected in the continuum limit. 

In the present analysis of the scaling of the maxima of the specific heat and suscep- 
tibility we found evidence for a cross-over regime in the four-dimensional Discrete 
Regge model of quantum gravity in the negative coupling region [23]. If this can be 
substantiated by further investigations and also for the Regge theory with contin- 
uous link lengths, the existence of a continuum limit at negative bare gravitational 
coupling is questionable. This is of major concern for a continuum theory of quan- 
tum gravity with matter fields [25]. 
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Figure 1: Expectation values of the average link lengths as a function of the gravi- 
tational coupling for (a) SRC, (b) GTM, and (c) Z 2 RM. 




Figure 2: Expectation values of the average curvature as a function of the gravita- 
tional coupling for (a) SRC, (b) GTM, and (c) Z 2 KM. All models exhibit a related 
phase structure and could lie in the same universality class. 
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Figure 3: Lattice average of the squared link length of type i as a function of the 
gravitational coupling for (a) SRC, (b) GTM, and (c) Z 2 RM. The symbol □ denotes 
the edges, o the face diagonals, A the body diagonals, and x the hyperbody diagonal 
of a triangulated hypercube. The dashed lines are to guide the eyes. 



0.01 



0.005 



cd 



-0.005 



-0.01 



0.2 



0.1 



□ p=-0.155 
o p=-0.157 
a p=-0.159 



1 



2 3 
d 



(a) 



(b) 



CD 







□ p=-0.110 
o (3= 0.1 1 1 
a p=-0.112 

-0.2 ' 1 



10 



CD 



-10 



1 



2 3 
d 



□ P=-4.668 



2 3 
d 



(c) 



Figure 4: Curvature correlation functions for (a) SRC, (b) GTM, and (c) Z 2 RM 
in the vicinity of f3~ . The curves show fits to the Yukawa ansatz (8) with fixed 
parameter a = 2. The resulting masses m are compiled in table 1. 




Figure 5: Histograms of the Regge action R t = A t St from simulations of the Z2- 
Regge Model. The distributions indicate a first-order transition around f3+ = 22.3 
(upper plot) and a continuous transition around f3~ = —4.665 (lower plot). 
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Figure 6: FSS of (a) the specific-heat maxima C max and (b) the susceptibility max- 
ima Xmax close to j3~ as a function of the lattice size L. In particular the behavior of 
Xmax is indicative for a cross-over rather than a true (continuous) phase transition. 




Figure 7: Volume dependence of the (pseudo-)critical gravitational couplings 
/3(C max ) and /3(x max ), as determined from the locations of C max and Xmax shown 
in figure 6. 
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